Looking for Light Pseudoscalar Bosons in the Binary Pulsar System J0737-3039 
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We present numerical calculations of the photon-light-pseudoscalar-boson conversion in the re- 
cently discovered binary pulsar system J0737-3039. Light pseudoscalar bosons (LPBs) oscillate into 
photons in the presence of strong magnetic fields. In the context of this binary pulsar system, 
this phenomenon attenuates the light beam emitted by one of the pulsars, when the light ray goes 
through the magnetosphere of the companion pulsar. We show that such an effect is observable in 
the gamma-ray band since the binary pulsar is seen almost edge-on, depending on the value of the 
LPB mass and on the strenght of its two-photon coupling. Our results are surprising in that they 
show a very sharp and significant (up to 50%) transition probability in the gamma-ray (> tens of 
MeV) domain. The observations can be performed by the upcoming NASA GLAST mission. 



Light pseudoscalar bosons (LPBs) coupled to two pho- 
tons have attracted considerable interest in the last few 
years and their implications for elementary-particle phe- 
nomenology, astrophysics and cosmology have been thor- 
oughly investigated 0, 0i 0- Specifically, the LPB 
mass m is assumed to be m < 1 eV and the LPB- 
photon-photon coupling is described by the interaction 
lagrangian 

A„* = -^f^-F,.0=I.E.B</>, (1) 

where </> denotes the LPB field, F'^'^ is the usual electro- 
magnetic field strength (F^^ is its dual), and M is the 
inverse coupling constant with the dimension of a mass. 
Natural Lorentz-Heaviside units with h = c ^ 1 are em- 
ployed throughout. 

Because the mass eigenstates of the LPB-photon sys- 
tem differ from the corresponding interaction eigenstates, 
interconversion takes place, a phenomenon quite similar 
in nature to flavour oscillations for massive neutrinos. 
However, since the off-diagonal term involves "two 
photons" , one of them actually corresponds to an exter- 
nal field. So, photon-LPB oscillations occur only in the 
presence of external electric or magnetic fields. 

A particular case of LPB is the standard axion, namely 
the pseudo-Goldstone boson associated with the Peccei- 
Quinn C/(l) global symmetry invented to solve the strong 
CP problem 0. Although the exact value of the ax- 
ion mass TO is model-dependent, generally one finds m ~ 
0.7 (10^° GeV/M) eV. Besides providing a natural candi- 
date for nonbaryonic dark matter the standard axion 
can in principle be detected with high-precision labora- 
tory experiments thanks to its two-photon coupling in 
eq. In this connection, two strategies have been pro- 



posed. One is based on the resonant axion-photon con- 
version inside a tunable microwave cavity jQ], while the 
other relies upon the induced ellipticity and change of 
polarization plane of an initially linearly-polarized laser 
beam JJ . In either case, a strong external magnetic field 
must be present. Generic LPBs with a two- photon cou- 
pling can also be discovered in this way, provided of 
course that their mass m and inverse couplig constant M 
fall into suitable ranges which depend on the experimen- 
tal arrangement. 

Quite recently, using an experimental setup based on 
ref. 01 , the PVLAS collaboration has observed an unex- 
plained effect that can be interpreted in terms of photon 
oscillations into a LPB with mass to ~ 1.0 • 10~^ eV and 
inverse coupling constant M ~ 3.8 • 10^ GeV 0|. Clearly 
these values prevent the LPB in question to be the stan- 
dard axion. 

On the other hand, LPBs are copiously produced in the 
core of stars, owing again to the two-photon coupling 
Therefore, under the assumption that LPBs escape from 
stars, energy-loss arguments put strong bounds on M. 
Typically, one gets M > IQi^GeV 0. This conclusion 
is further supported by the result of the CAST experi- 
ment Manifestly, the PVLAS result contradicts this 
bound by about five orders of magnitude. A possible way 
out of this difficulty is to suppose that, because of some 
new physics, the LPB detected by PVLAS is trapped in- 
side star cores. A specific scenario in which this indeed 
occurs has been put forward 0]. While no fully satis- 
factory solution to this problem has yet been found, the 
existence of new physics at low energy is evidently called 
for, and this poses an exciting challenge for contempo- 
rary research. Thus, the need for an independent test of 
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the PVLAS result looks imperative. 

In this Letter, we show that LPBs can give rise to 
observable effects if suitably selected binary pulsars are 
monitored in the gamma-ray band with future satellite- 
borne detectors like GLAST. Our analysis is carried out 
for arbitrary values of m and M, and the requirement of 
observability defines the corresponding parameter ranges. 
It turns out that the values found by PVLAS fall within 
these ranges, and so our result provides a cross-check for 
the PVLAS claim of a new LPB. Actually, our calcula- 
tions apply also to the case of scalar bosons, except that 
the E • B coupling should be replaced by the E ■ E — B • B 
coupling. 

In what follows we will focus, as a model system, on 
the specific case of the recently-discovered binary pul- 
sar system J0737-3039 which is seen almost edge on ^3 
and for which gravitational microlensing has also been 
studied We have already shown .13] that in the pe- 
culiar geometry of this system it will be possible to ob- 
serve quantum vacuum lensing resulting from the photon- 
magnetic field interaction, with photon energies in the 
X-ray domain. There we used the same system, again 
starting with the photon beam from one of two stars 
(pulsar A). Denote by p the impact parameter of the 
beam, namely its projected distance from the compan- 
ion, pulsar B. When p is small, the beam goes through 
the magnetosphere of pulsar B. Therefore the observed 
photons traverse a region where a strong magnetic field 
is present. Assume now that a LPB exists. Then, owing 
to its two-photon coupHng (1), some of the photons will 
be converted into LPBs, thereby decreasing the beam in- 
tensity. Depending on the actual values of the physical 
parameters involved, the beam attenuation can be sub- 
stantial. As we will see, the effect turns out to be impor- 
tant for photons with energy larger that 10 MeV, namely 
in the gamma-ray band. Remarkably enough, these pho- 
tons propagate totally unimpeded in the magnetosphere 
of pulsar B, and so we do not have to bother about fur- 
ther potential sources of beam attenuation. Moreover, 
the effect under consideration does not depend on any 
new physics responsible for the confinement of the LPBs 
in the star interiors . 

The model system studied, J0737 — 3039, is particu- 
larly well suited for our purposes. Its orbital plane has 
currently an inclination angle i ~ 87° ^■'^d the or- 
bital precession of the system allows us to predict that 
the value i ~ 90" could be achieved before 2020. The 
orbital period is 2 hours and 45 minutes, while the spin- 
ning period of pulsar B is 2.77 seconds. Although the 
current minimum value of the impact parameter of the 
light beam from pulsar A is somewhat uncertain, a re- 
alistic estimate yields p ~ 4 • 10'^ km and this is the 
value used in our analysis. The geometry of the binary 
pulsar under consideration is depicted in Fig. 

The light beam from pulsar A propagates along the 
z direction and passes pulsar B at impact parameter p. 
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FIG. 1: Geometry of the model binary pulsar system 

We introduce a (x, y, z) coordinate system fixed with pul- 
sar B. For simplicity, we suppose that the rotation axis 
of pulsar B lies along the y direction and that its mag- 
netic dipolc moment precedes at an angle 9 = 45" about 
the y axis. The light beam from pulsar A experiences 
a time-dependent and inhomogcncous dipolar magnetic 
field B. The photon-LPB oscillations are described by 
the coupled Klein- Gordon and Maxwell equations dic- 
tated by lagrangian (^), plus free-field terms. In the 
short- wavelength approximation (to be justified later), 
the propagation equations take the following first-order 
form [3 

, = Mmz)) , (2) 

oz 

with 

\tp(z)) = c^{z) \x) + Cy{z) \y) + C0 , (3) 

where \x) and \y) are the two linear polarization states of 
the photon (propagating in the z-direction) along the x 
and y axis, respectively, and \(j)) is the LPB state. More- 
over, the mixing matrix M. is given by 

fuj + A,, A,,y BJ2M \ 
yx ^ ~^ ^yy 

By/2M , (4) 
\ B^/2M By/2M LO - m^/2uj J 

with uj denoting the photon/LPB energy. While the 
terms appearing in the third row and column of Ai have 
an evident physical meaning, the A-terms require fur- 
ther discussion. Generally speaking, they reflect both the 
properties of the medium and the QED vacuum polar- 
ization effects. Off-diagonal terms A.xy, Ayx directly mix 
the photon polarization states and typically give rise to 
Faraday rotation. Clearly they are important when the 
polarization of photon beam is taken into account, but 
become irrelevant for the present analysis. The diago- 
nal elements A^x , Ayy arise from two independent contri- 
butions: vacuum polarization effects and plasma effects. 
The first ones are due to one-loop Feynman diagrams and 
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are analytically described by the Heisenberg-Euler effec- 
tive lagrangian when terms involving LPBs in the loops 
are discarded, as is currently done. Neglecting again 
beam polarization effects, we get A^^f^ ~ ^yy"^ — 
{a / ibTr) {B / Bcr)"^ uj, where a is the fine-structure con- 
stant and Br-r ~ 4.41 • lO^'^ G denotes the critical mag- 
netic field [l5j . Plasma effects arise from the presence of 
free charges in the magnetosphere of pulsar B and their 
contribution is A^'^ ~ A£!f ~ — 27ra n„ luj m,„ , where 
is the electron mass and rig is the electron density |14| . 
Off-diagonal QED effects and plasma effects turn out to 
be presently irrelevant and will be neglected (more about 
this, later). 

We solve eq. Q numerically in the equivalent expo- 
nential form 



\i^{z)) =Tcxpi-i M{z')dz' l^(zo)) 



(5) 



where T is Dyson's ordering operator. Since the ma- 
trix is explicitly z-dependent, the propagation is com- 
puted iteratively, i.e. in steps of small intervals 5z. For 
small enough 5z, the M. matrix can be considered as con- 
stant over one distance step, and the short-propagation 
iterations can be calculated as 



\ij}{z + 5z)) ~ exp {-iM{z)5z)\^:{z)) . 



(6) 



The propagation is performed from zq to Zmax with the 
initial condition corresponding to a pure photon state. 
Here zq is fixed by the trajectory of pulsar A and Zmax is 
chosen so that Bx{z)/2M ~ By{z)/2M ~ for z > z^ax- 
The transition probability is then 



= \{4>\i'{Zmax)}f = \c4,{Zjnax)\^ 



(7) 



where j = x (resp. y) refers to the initial photon polariza- 
tion along the x-axis (resp. y-axis). For an unpolarized 
photon beam, the transition probability is 



Pi 



(7^0) 



P, 



(8) 



In what follows we use the values of m and M as 
suggested by the PVLAS result. In Fig. 2, we have 
plotted the transition probability P(^^0-) versus the pho- 
ton energy w, for a gamma beam with impact parameter 
p = 4 • 10^ km. 

The strong energy-dependence exhibited by P(^^0) 
suggests to investigate the photon-LPB conversion in the 
region of the broad maximum, namely for lu > lOMeV. 
Observe that this circumstance fully justifies the above 
short-wavelength approximation. 

Actually, the threshold behaviour of Pj..^^^^) can be 
understood as follows. Consider a drastically simpli- 
fied situation in which the magnetic field experienced 
by the light beam is stationary and homogeneous. A 
typical value of the magnetic field strenght on the sur- 
face of a neutron star, namely at a distance r ~ 10 km 
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FIG. 2: Transition probability versus photon energy for a 
trajectory of the light beam with an impact parameter p = 
4 • 10^ km. 



from the centre, is Bq ~ 10^^ G. Since B is dipolar, it 



goes like 



and so the magnetic field strenght at 



p ~ 4 ■ 10'^ km is P ~ 1.6 • 10^ G. Correspondingly, we 
have AQ^^ ~ A^/^^ ~ 6.3- IQ-i® (uj/MeV) eV. Because 
the terms Aj.y, Aj,^, presently get only a QED contribu- 
tion, they should not exceed A^^-^ ~ A'^^^. Follow- 
ing ref. [iJI, it is straightforward to see that the terms 
AP^ ~ A^jf are much smaller than A^.^^ ^ A^^f ^ 
for P ~ 1.6 • lO^G. We find P^/2A/ ~ Pj,/2M ~ 
0.4 • 10-12 eV and m^/2uj ~ 0.5 • IQ-^^ (MeV/w) eV. Cor- 
respondingly Axy and Ayx turn out to be the smallest el- 
ements oi M. Thus, we see that off-diagonal QED effects 
as well as plasma effects are indeed negligible (in agree- 
ment with previous assumptions). In addition, mixing 
effects in A4 should become physically important when 
the ratio of the off-diagonal elements to the difference 
of the diagonal ones starts to be of order 1. Therefore, 
photon-LPB conversion is expected to actually occur for 
LU > 1 MeV, which means that a substantial beam atten- 
uation ought to show up when gamma-ray photons above 
10 MeV are observed. 

It makes even more sense, at this point, to consider 
the case of J0737, because pulsars are frequently celes- 
tial gamma-ray sources ^iJl , because of the very special 
geometry of this system and because the NASA GLAST 
mission (due to launch in 2007) will have a pulsar gamma- 
ray detection sensitivity up to two orders of magnitude 
better than any previous telescope [l^ . 

Assuming that the gamma photon beam emitted by 
pulsar A is not polarized, the temporal behaviour of the 
considered effect is best expressed in terms of the total 
transmission A — 1 — Pj.^^^) of the gamma ray after 
propagation in the magnetosphere of pulsar P. In Fig. 
3, we plot A versus time as pulsar A moves in its nearly 
edge-on orbit. Our numerical simulation predicts in this 
case a strong attenuation of the photon beam up to 50 % 
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FIG. 3: Total transmisson of the gamma photon beam emit- 
ted by pulsar A versus time. Inset shows the modulation 
mainly due to the rotation of the magnetic dipole moment of 
pulsar B. 



with a time duration of about 200 s. The effect has three 
different temporal structures. Of course, the broad min- 
imum, from —200 s to -f 200 s, corresponds to the transit 
of pulsar B "in front" of pulsar A. The tens-of-seconds, 
symmetric peaks are due to photon-LPB oscillations de- 
pending on the path through the interaction region with 
pulsar B. Finally, the highest frequency modulation, 
shown in inset, is due to the rotation of the magnetic 
dipole moment of pulsar B. Under the assumptions cor- 
responding to the PVLAS result, a macroscopic effect 
could be observable from the J0737-3039 system, given 
enough instrumental sensibility and observing time. We 
note that pulsar A also generates a similar effect on the 
gamma photon beam emitted by pulsar B. Since the 
spin period of pulsar A is smaller than the one of pulsar 
B (23 ms), the modulation period of the effect will be 
also smaller. 
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FIG. 4; Exclusion region in the case that the existence of the 
attenuation is excluded at 10 % level. 



Finally, Fig. 4 shows the region in the m — M pa- 
rameter space which is excluded by no detection of an 
attenuation A at the 10 % level. To get such an atten- 
uation, one need to count 100 photons during the total 
integrating time. For two weeks observation time, this 
corresponds to a flux from pulsar A of about 2 • 10~^ 
photons/cm^ /s which is a reasonable flux according to 
previous observations of several pulsars • We note in 
any case that the GLAST sensibility curves allow for a 
much weaker minimum detectable flux, down to 1 • 10"^" 
photons/cm^/s In conclusion, we have shown that 
LPBs can produce observable effects in a binary pulsar 
system seen almost edge-on, like J0737-3039. The effect 
is in the gamma-ray band, depending of course on the 
actual values of the parameters m and M . This is an in- 
triguing feature of our result, since pulsars are frequently 
celestial gamma-ray sources and the NASA GLAST mis- 
sion (due to launch in 2007) will have a pulsar gamma-ray 
detection sensitivity up to two orders of magnitude better 
than any previous telescope |l8l |. Of course, new binary 
pulsar systems will be discovered soon through radio- 
astronomy or possibly by the very GLAST mission. A 
potentially useful new tool seems to be coming of age. 
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